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Total absorption spectroscopy was used to investigate the /?-decay intensity to states above the 
neutron separation energy followed by y-ray emission in and ®^Rb. Accurate results were 

obtained thanks to a careful control of systematic errors. An unexpectedly large 7 intensity was 
observed in all three cases extending well beyond the excitation energy region where neutron pen¬ 
etration is hindered by low neutron energy. The 7 branching as a function of excitation energy 
was compared to Hauser-Feshbach model calculations. For ®^Br and ®®Br the 7 branching reaches 
57% and 20% respectively, and could be explained as a nuclear structure effect. Some of the states 
populated in the daughter can only decay through the emission of a large orbital angular momen¬ 
tum neutron with a strongly reduced barrier penetrability. In the case of neutron-rich ®‘*Rb the 
observed 4.5% branching is much larger than the calculations performed with standard nuclear sta¬ 
tistical model parameters, even after proper correction for fluctuation effects on individual transition 
widths. The difference can be reconciled introducing an enhancement of one order-of-magnitude in 
the photon strength to neutron strength ratio. An increase in the photon strength function of such 
magnitude for very neutron-rich nuclei, if it proved to be correct, leads to a similar increase in the 
(n,7) cross section that would have an impact on r process abundance calculations. 

PACS numbers: 23.40.-s, 21.10.PC, 29.30.Kv, 26.50.-l-x 


Neutron unbound states can be populated in the /3 de¬ 
cay of very neutron-rich nuclei, when the neutron separa¬ 
tion energy S'„ in the daughter nucleus is lower than the 
decay energy window Qfj. Given the relative strengths 
of strong and electromagnetic interactions these states 
decay preferentially by neutron emission. Beta delayed 
7-ray emission from states above Sn was first observed 
in 1972 in the decay of ®^Br a. Since then it has been 
observed in a handful of cases: 0, ®^Rb 0, ®®As Q, 

141 Cs a, a, and 94Rb 0. The paucity of infor¬ 

mation is related to the difficulty of detecting weak high- 
energy 7-ray cascades with the germanium detectors that 
are usually employed in /-decay studies. This problem 
has become known as the Pandemonium effect it 

also affects the accuracy of the data. 

There is an analogy Q between this decay process and 
neutron capture reactions which populate states in the 
compound nucleus that re-emit a neutron (elastic chan¬ 
nel) or de-excite by 7 rays (radiative capture). Indeed the 
reaction cross section is parametrized in terms of neutron 
and 7 widths, r„ and E..^ respectively, which also deter¬ 


mines the fraction of / intensity above Sn that proceeds 
by neutron or 7 emission. Radiative capture (n, 7) cross 
sections for very neutron-rich nuclei are a key ingredient 
in reaction network calculations used to obtain the yield 
of elements heavier than iron in the rapid (r) neutron 
capture process occurring in explosive-like stellar events. 
It has been shown |inl4l2| that the abundance distribu¬ 
tions in different astrophysical scenarios are sensitive to 
(n, 7) cross sections. In the classical “hot” r process late 
captures during freeze-out modify the final element abun¬ 
dance. In the “cold” r process the competition between 
neutron captures and / decays determines the forma¬ 
tion path. Cross section values for these exotic nuclei 
are taken from Hauser-Feshbach model calculations [l^ , 
which are based on a few quantities describing average 
nuclear properties: nuclear level densities (NLD), pho¬ 
ton strength functions (PSF) and neutron transmission 
coefficients (NTC). Since these quantities are adjusted to 
experiment close to j3 stability it is crucial to find means 
to verify the predictions for very neutron-rich nuclei. 

The Total Absorption Gamma-ray Spectroscopy 
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(TAGS) technique aims at detecting cascades rather than 
individual 7 rays using large 47 r scintillation detectors. 
The superiority of this method over high-resolution ger¬ 
manium spectroscopy to locate missing (3 intensity has 
been demonstrated before 14, El- However its appli¬ 
cation in the present case is very challenging, since the 
expected 7 -branching is very small and located at rather 
high excitation en erg ies. As a matter of fact previous 
attempts at LNPI [7| with a similar aim did not lead to 
clear conclusions. In this Letter we propose and demon¬ 
strate for the first time the use of the TAGS technique 
to study 7 -ray emission above S'„ in /3-delayed neutron 
emitters and extract accurate information that can be 
used to improve (n, 7 ) cross section estimates far from (3 
stability. 

Neutron capture and transmission reactions have been 
extensively used [l^ to determine neutron and 7 widths 
(or related strength functions). An inspection of Ref. 161 
shows that in general r„ is orders-of-magnitude larger 
than r.^. In the decay of ®’^Br, which is the best stud¬ 
ied case a dozen states emitting single 7 rays 

were identified within 250 keV above Sn collecting about 
0.5% of the decay intensity to be compared with a neu¬ 
tron emission probability of 2.6%. The observation of 
such relatively high 7 -ray intensity was explained as be¬ 
ing due to a nuclear structure effect: some of the levels 
populated can only decay by emission of neutrons with 
large orbital angular momentum I, which is strongly hin¬ 
dered. In addition it has been pointed out that a siz¬ 
able 7 -ray emission from neutron unbound states can be 
a manifestation of Porter-Thomas (PT) statistical fluc¬ 
tuations in the strength of individual transitions. The 
role and relative importance of both mechanisms should 
be investigated. 

We present here the results of measurements for 
three known neutron emitters, ®^Br 2l|, ®®Br 2^ and 
94Rb [13 , using a newly developed TAGS spectrome¬ 
ter. The results for ®^Rb, also measured, will be pre¬ 
sented later 2^. The measurements were performed 


at the IGISOL mass separator [23 of the University 
of Jyvaskyla. The isotopes were produced by proton- 
induced fission of uranium and the mass-separated 
beam was cleaned from isobaric contamination using the 
JYFLTRAP Penning trap 2^ 271. The resulting beam 
was implanted at the centre of the spectrometer onto a 
movable tape which periodically removed the activity to 
minimize daughter contamination. Behind the tape was 
placed a 0.5 mm thick Si detector with a /3-detection 
efficiency of about 30%. The Valencia-Surrey Total Ab¬ 
sorption Spectrometer Rocinante is a cylindrical 12-fold 
segmented BaF 2 detector with a length and external di¬ 
ameter of 25 cm, and a longitudinal hole of 5 cm diam¬ 
eter. The detection efficiency for single 7 rays is larger 
than 80%. The spectrometer has a reduced neutron sen¬ 
sitivity in comparison to Nal(Tl) detectors, a key feature 
in the present application. It also allows the measure¬ 


ment of multiplicities which helps in the data analysis. 
In order to eliminate the detector intrinsic background 
and the ambient background we use /3-gated TAGS spec¬ 
tra in the present analysis. Nevertheless other sources of 
spectrum contamination need to be characterized accu¬ 
rately. 

In the first place the decay descendant contamination, 
was computed using the Geant4 simulation toolkit ( 2 ^ . 
In the case of the daughter decay we use an event gener- 
ator based on the well known decay level scheme 
The calculated normalization factor was adjusted to pro¬ 
vide the best fit to the recorded spectrum. The measure¬ 
ment of ®®Br was accidentally contaminated by ®'^Y, the 
long-lived grand-daughter of ®^Rb, and was treated in 
the same manner. The case of the contamination due to 
the /3-delayed neutron branch is more challenging. The 
decay simulation must explicitly include the emitted neu¬ 
trons. These neutrons interact with detector materials 
producing 7 rays through inelastic and capture processes. 
An event generator was implemented which reproduces 
the known neutron energy distribution, taken from [ 2 ^ . 
and the known 7 -ray intensity in the final nucleus, taken 
from 21-^. The event generator requires the /3 inten¬ 
sity distribution followed by neutron emission which 
was obtained from deconvolution of the neutron spec¬ 
trum. Another issue is whether the interaction of neu¬ 
trons with the detector can be simulated accurately. We 
have shown recently [s^ that this is indeed the case pro¬ 
vided that Geant4 is updated with the newest neutron 
data libraries and the original capture cascade generator 
is substituted by an improved one. The normalization 
factor of the /3-delayed neutron decay contamination is 
fixed by the Pn value. Another important source of spec¬ 
trum distortion is the summing-pileup of events. If more 
than one event arrives within the same ADG event gate, 
a signal with the wrong energy is stored in the spectrum. 
Apart from the electronic pulse pile-up effect for a single 
detector module [3l| one must consider the summing of 
signals from different detector modules. A new Monte 
Carlo (MG) procedure to calculate their combined con¬ 
tribution has been developed. The procedure is based 
on the random superposition of two stored events within 
the ADC gate length. The normalization of the resulting 
summing-pileup spectrum is fixed by the event rate and 
the ADC gate length 311. 

Several laboratory 7 -ray sources were used to deter¬ 
mine the energy and resolution calibration of the spec¬ 
trometer. The measured singles spectra also served to 
verify the accuracy of the spectrometer response simu¬ 
lated with Geant4. The use of /3-gated spectra in the 
analysis required additional verifications of the simula¬ 
tion. Due to the existence of an electronic threshold in 
the Si detector (100 keV) the /3-detection efficiency has a 
strong dependence with /3-endpoint energy up to about 
2 MeV. This affects the region of interest (see Fig.[T]). To 
verify that the MG simulation reproduces this energy de- 
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FIG. 1. Relevant histograms for ®®Br: parent decay (grey 
filled), daughter decay (pink), summing-pileup (green), f}- 
delayed neutron decay (dark blue), accidental contamination 
(light blue), reconstructed spectrum (red). 


pendence we use the information from a separate experi¬ 
ment 32| measuring P„ values with the neutron counter 


BELEN and the same /3 detector. Several isotopes with 
different neutron emission windows Qp — Sn were mea¬ 
sured, resulting in variations of the neutron-gated /3 ef¬ 
ficiency as large as 25%. Geant4 simulations using the 
above mentioned /3-delayed neutron decay generator are 
able to reproduce the isotope-dependent efficiency within 
better than 4%. 


Figure [U shows the /3-gated TAGS spectrum measured 
during the implantation of ®®Br ions. Also shown is the 
contribution of the daughter ®®Kr decay, the neutron de¬ 
cay branch populating ®^Kr, the summing-pileup contri¬ 
bution and the accidental contamination of Notice 
the presence of net counts beyond the neutron separation 
energy, which can only be attributed to the decay feeding 
excited states above Sn which de-excite by y-ray emis¬ 
sion. In this region the major background contribution 
comes from summing-pileup which is well reproduced by 
the calculation as can be observed. Similar pictures were 
obtained for the decay of ®^Br and ®^Rb. 

The analysis of the /3-gated spectra follows the method 
developed by the Valencia group [H, in. The inten¬ 
sity distribution is obtained by deconvolution of the 
TAGS spectrum with the calculated spectrometer re¬ 
sponse to the decay. The response to electromagnetic 
cascades is calculated from a set of branching ratios (BR) 
and the MG calculated response to individual 7 rays. 
Branching ratios are taken from (ill - E^ for the low en¬ 
ergy part of the decay level scheme. The excitation en¬ 
ergy range above the last discrete level is treated as a con¬ 
tinuum divided into 40 keV bins. Average BR for each 
bin are calculated from NLD and PSF as prescribed by 
the Hauser-Feshbach model. We use NLD from Ref. [35| 
as tabulated in the RIPL-3 library [s^. The PSF is ob¬ 
tained from Generalized Lorentzian (El) or Lorentzian 


(Ml, E2) functions using the parameters recommended 
in Ref. [Sg. The electromagnetic response is then con¬ 
voluted with the simulated response to the /3 continuum. 
The spin-parity of some of the discrete states at low exci¬ 
tation energy in the daughter nucleus is uncertain. They 
are however required to calculate the BR from the states 
in the continuum. The unknown spin-parities were varied 
and those values giving the best reproduction of the spec¬ 
trum were adopted. There is also ambiguity in the spin- 
parity of the parent nucleus which determines the spin- 
parity of the levels populated in the continuum. Here we 
assume that allowed Gamow-Teller (GT) selection rules 
apply. Our choices, 3/2“ for ®^Br, 1“ for ®®Br and 3“ 
for ®'^Rb, are also based on which values best reproduce 
the spectrum. 



(red), high-resolution 7 spectroscopy (blue), from /3-delayed 
neutron (grey filled). 


As an example of the results of the analysis we show 
in Fig. [2] the Ip-y intensity obtained for ®®Br. The spec¬ 
trum reconstructed with this intensity distribution re¬ 
produces well the measured spectrum (see Fig. [T|). The 
analysis for the other two isotopes shows similar quality 
in the reproduction of the spectra. We also include in 
Fig. [5] the intensity obtained from high-resolution mea¬ 
surements [ 2 ^, showing a strong Pandemonium effect. 
The Pandemonium effect is even stronger in the case of 
®'*Rb and somewhat less for ®^Br. The complete and 
its impact on reactor decay heat 37] and antineutrino 
spectrum [33 summation calculations will be discussed 
elsewhere [39] . Here we concentrate on the portion of 
that intensity located in the neutron unbound region. A 
sizable TAGS intensity is observed above Sn extending 
well beyond the first few hundred keV where the low neu¬ 
tron penetrability makes 7 -ray emission competitive. For 
comparison Fig.J^also shows Ipn deduced from the neu¬ 
tron spectrum [29] as explained above. The Ipj above 
Sn adds up to = 1.59(17)%, to be compared with 

the integrated Ipn (or P„) of 6.4(6)%. From the TAGS 
analysis for the other two isotopes we find a ^ Ip-y of 
3.5(3)% (®^Br) and 0.53(14)% (®^Rb) to be compared 
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with P„-values of 2.60(4)% and 10.18(24)% respectively. 
In the case of ®^Br we find 7 times more intensity than 
the high-resolution measurement [l9l |. The uncertainty 
quoted on Ip-y is dominated by systematic uncertain¬ 
ties. We did a careful evaluation of possible sources of 
systematic effects for each isotope. The uncertainty com¬ 
ing from assumptions in the BR varies from 1% to 5% 
(relative value) depending on the isotope. The impact 
of the use of different deconvolution algorithms is in 
the range of 2% to 10%. The uncertainty in the energy 
dependence of the /3 efficiency contributes with 4%. The 
major source of uncertainty comes from the normaliza¬ 
tion of the background contribution, which at the ener¬ 
gies of interest is dominated by the summing-pileup. We 
estimated that reproduction of spectra could accommo¬ 
date at most a ±15% variation from the nominal value, 
which translates into uncertainties of 6 % to 22 %. 

Figure [3] shows the ratio Ip-y/{Ip-y +1pn) in the range of 
energies analyzed with TAGS for all three cases. This ra¬ 
tio is identical to the average ratio (r.y/(r.y±r„)) over all 
levels populated in the decay. The shaded area around 
the experimental value in Fig. [3] serves to indicate the 
sensitivity of the TAGS results to background normal¬ 
ization as indicated above. The average width ratio was 
calculated using the Hauser-Feshbach model. The results 
for the three spin-parity groups populated in GT decay 
are shown. The NLD and PSF values used in these cal¬ 
culations are the same as those used in the TAGS anal¬ 
ysis. The new ingredient needed is the NTG, which is 
obtained from the Optical Model (OM) with the TALYS- 
1.4 software package [i^. OM parameters are taken from 
the so-called local parametrization of Ref. ^ . Neutron 
transmission is calculated for known final levels popu¬ 
lated in the decay 21-2^. In order to compute the av¬ 


erage width ratio we need to include the effect of statis¬ 
tical fluctuations in the individual widths [i^. We use 
the MG method to obtain the average of width ratios. 
The sampling procedure is analogous to that described 
in Ref. (33l| . Level energies for each spin-parity are gener¬ 
ated according to a Wigner distribution and their corre¬ 
sponding F..^ and r„ to individual final states are sampled 
from PT distributions. The total 7 and neutron widths 
are obtained by summation over all possible final states 
and the ratio computed. The ratio is averaged for all 
levels lying within each energy bin. In order to suppress 
fluctuations in the calculated average, the sampling pro¬ 
cedure is repeated between 5 and 1000 times depending 
on level density. Very large average enhancement factors 
were obtained, reaching two orders-of-magnitude when 
the neutron emission is dominated by the transition to a 
single final state. 

In the case of ®^Br 3/2“ decay one can see in Fig. [3] 
that the strong 7 -ray emission above Sn can be explained 
as a consequence of the large hindrance of Z = 3 neutron 
emission from 5/2“ states in ®^Kr to the O’*" g.s. of ®®Kr, 
as pointed out in Ref. [i[. In the case of ®®Br 1“ decay a 
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FIG. 3. Average gamma to total width ratio from experiment 
and calculated for the three spin-parity groups populated in 
allowed decays. The shaded area around the experimental 
value indicates the sensitivity to the background normaliza¬ 
tion (see text). 


similar situation occurs for 0“ states in ®®Kr below the 
first excited state in ®^Kr at 532 keV, which require I = 3 
to populate the 5/2+ g.s. in ®’^Kr. For a more quanti¬ 
tative assessment one should know the distribution of (3 
intensity between the three spin groups, which could be 
obtained from /3-strength theoretical calculations. The 
case of ^'^Rb 3“ decay is the most interesting. The fi¬ 
nal nucleus ®^Sr is five neutrons away from /3 stability. 
The 7 intensity although strongly reduced, only 5% of 
the neutron intensity, is detectable up to 1.5 MeV be¬ 
yond Sn- The structure observed in the average width 
ratio, is associated with the opening of /3n channels to 
different excited states. Note that the structure is re¬ 
produced by the calculation, which confirms the energy 
calibration at high excitation energies. In any case the 
calculated average gamma-to-total ratio is well below the 
experiment. In order to bring the calculation to the ex¬ 
perimental value one would need to enhance the PSF, 
or suppress the NTG, or any suitable combination of the 
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two, by a very large factor. For instance we verified that 
a twenty-fold increase of the El PSF would reproduce 
the measurement assuming a /3-intensity spin distribu¬ 
tion proportional to 2J -\- 1. An enhancement of such 
magnitude for neutron-rich nuclei, leading to a similar 
enhancement of (n, 7 ) cross sections, will likely have an 
impact on r-process abundance calculations. Therefore 
it will be important to investigate the magnitude of pos¬ 
sible variations of the NTC. 

In conclusion, we have confirmed the suitability of the 
TAGS technique to obtain accurate information on 7 -ray 
emission from neutron unbound states and applied it to 
three known /3-delayed neutron emitters. A surprisingly 
large 7 -ray branching of 57% and 20% was observed for 
®^Br and ®®Br respectively, which can be explained as a 
nuclear structure effect. In the case of ®^Br we observe 7 
times more intensity than previously detected with high 
resolution 7 -ray spectroscopy, which confirms the need of 
the TAGS technique for such studies. In the case of the 
more neutron-rich ®^Rb the measured branching is only 
4.5% but still much larger than the results of Hauser- 
Feshbach statistical calculations, after proper correction 
for individual width fluctuations. The large difference 
between experiment and calculation can be reconciled by 
an enhancement of standard PSF of over one order-of- 
magnitude. To draw more general conclusions it will be 
necessary to extend this type of study to other neutron- 
rich /3-delayed neutron emitters. Such measurements us¬ 
ing the TAGS technique are already underway and addi¬ 
tional ones are planned. 
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